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CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] This appKcation claims priority, tinder 35 U.S.C. §119(e), of provisional 

application No. 60/280,793, filed April 2, 2001. 



STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR 
DEVELOPMENT 



[0002] Not applicable. 



BACKGROUND OF THE INVENTION 



[0003] This invention is in the field of digital communications, and is more 

specificaEy directed to Digital Subscriber Line (DSL) commimications. 



[0004] An important and now popular modulation standard for DSL 

5 communication is Discrete Multitone (DMT). According to DMT technology, the 
available spectrum is subdivided into many subchannels (e.g., 256 subchannels of 4.3125 
kHz). Each subchannel is centered about a carrier frequency that is phase and 
amplitude modulated, typically by Quadrature AmpHtude Modulation (QAM), in which 
each symbol value is represented by a point in the complex plane; the number of 
10 available symbol values depends, of course, on the number of bits in each sj^mbol. 
During initialization of a DMT communications session, the number of bits per symbol 
for each subchaimel (ie., the "bit loading") is determined according to the noise 
currently present in the transmission channel at each subchannel frequency and 
according to the transmit signal attenuation at that frequency. For example, relatively 
15 noise-free subchannels may communicate data in ten-bit to fifteen-bit symbols 
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corresponding to a relatively dense QAM constellation (with short distances between 
points in the constellation), while noisy channels may be limited to only two or three bits 
per symbol (to allow a greater distance between adjacent points in the QAM 
constellation). Indeed, some subchannels may not be loaded with any bits, because of 
5 the noise and attenuation in those channels. In this way, DMT maximizes the data rate 
for each subchannel for a given noise condition, permitting high speed access to be 
carried out even over relatively noisy twisted-pair lines. 

[0005] DMT modxilation also permits much of the processing of the data to be 

carried out in the digital domain. Typically, the incoming bitstream is seriaUy received 
^" 10 and then arranged into symbols, one for each subchannel (depending on the bit 
5 loading). Reed-Solomon coding and otiier coding techniques are also typically appHed 

^ for error detection and correction. Modulation of the subchannel carriers is obtained by 

ffl|5Sl 

# application of an inverse Discrete Fourier Transform (IDFT) to the encoded symbols, 

H producing the output modulated time domain signal. This modulated signal is then 

15 serially transmitted. All of these operations in DMT modulation can be carried out in 
the digital domain, permitting implementation of much of a DSL modem, and 
particularly much of the processing-intensive operations, in a single chip (such as a 
Digital Signal Processor, or DSP). 

[0006] The discrete output time domain signal from the modulation is Hciea 

20 converted into a time-domain analog signal by a conventional digital-to-analog 
converter. The analog signal is then communicated over the transmission channel to the 
receiving modem, which reverses the process to recover the transmitted data. The non- 
ideal impulse response of the trai«mission chaimd of course distorts the transmitted 
signal. Accordingly, the signal received by the receiving modem wiU be a convolution 
25 of the analog output waveform with the impulse response of the transmission channel. 
Ideally, the DMT subchannels in the received signal are orthogonal so that the 
modulating data can be retrieved from the transmitted signal by a Discrete Fourier 
Transform (DFT) demodtdation, under the assumption that convolution in the time 
domain corresponds to multiplication in the frequency domain. 
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[0007] While DMT provides excellent transmission data rates over modest 

comtnuiucations facilities such as twisted-pair wires, the IDFT modulation can resxilt in 
a high peak-to-average ratio (PAR) of the signal amplitudes. The PAR is defined as the 
ratio of the peak power level, for a sample, to the average power level over a sequence of 

5 samples. In the case of conventional DMT modulation, the amplitude of the time 
domain signal from the IDFT has a probability distribution function that has 
substantially a Gaussian shape. This Gaussian distribution of the signal amplitudes 
indicates the possibility that some time-domain samples may have amplitudes that are 
very high, as compared with the average sample amplitude. The resulting PAR is 

10 therefore much higher for DMT signals than for single-channel signals, because of the 
probability that the additive DMT peaks can overlay one another to result in an 
extremely large amplitude time-domain sample. 



;p [0008] The high PAR for conventional DMT signals presents significant 

W constraints on the transmission circuitry, and can greatly compUcate the analog circuitry 

r 15 required for high fidelity transmission. For example, a high PAR translates into a large 
^ dynamic range at the inputs of digital-to-analog and analog-to-digital converters, 

necessitating a large number of bits of resolution, and the associated extreme cost and 
complexity. Filters and amplifiers must also become more complex and costly in order 
to handle both the high peak ampUtudes and also the resolution required for the vast 
20 majority of the samples having lower ampUhide. hi addition, the high PAR results in 
much higher power consumption in the communications circuits, further increasing the 
cost of the circuits and systems used for DMT transmission and receipt, particularly 
those circuits often referred to as the analog front end (AFE)v 

[0009] A common approach to controlling the PAR in DMT transmission is to 

25 cHp amphfaides that exceed a selected threshold. This clipping obviously results in loss 
of signal. Clipping effectively introduces an impulse at the clipped sample in the time 
domain signal, having the negative of the ampHtude being cUpped. As known in the art, 
a time-domain impulse corresponds to additive noise across aU subchannels in the 
frequency domain, and thus clipping effectively reduces the signal-to-noise ratio for all 
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subchannels in the modulated signal, hi addition, the use of cHpping requires a difficult 
tradeoff. Clipping to a relatively low ampUtude threshold reduces circuit complexity 
and power dissipation, but greatly increases the probability of dipping; on the other 
hand, clipping at a high amplitude threshold decreases the probability of cHpping but 
5 increases circuit complexity and cost. If the average power of the signal is kept small, so 
that the peak ampUtude remains within the d3mainic range of the analog circuitry, the 
signal amplifiers are operating in an inefficient operating state; conversely, the amplifier 
efficiency can be improved by raising the average power only by increasing the 
probability of clipping. 

10 [0010] Various approaches have thus been developed to reduce the PAR of DMT 

O signals to minimize the number of samples that require clipping. According to one class 

y of techniques, DMT symbols are coded so that the resulting code words reside in the set 

% of DMT symbols that reside below the desired PAR amplitude threshold. These 

U1 techniques necessitate a loss of data rate, because of the coding overhead that results. 

L 15 [0011] Another approach effects an invertible transformation on the transmitted 

|! sigrial, such as a phase rotation for certain subchannels, to reduce the probabiUty that the 

PAR amplitude threshold is exceeded. Assuming the probabiHty of the original signal 
exceeding tiie PAR threshold to be low, the probabihty that both the original signal and 
also the traiwformed signal will exceed the threshold will be approximately the square 
20 of the low probabiHty for the original signal, which greatly reduces the extent of 
cUpping. hi this approach, a control signal is commtmicated to the receiver to identify 
the transformation, so that the receiver can apply the inverting transformation, as 
necessary, and recover the original signal 

[0012] Another approach estimates and corrects the effects of clipping at the 

25 receiver. As described in European Patent Application publication EP 0957 615 A2, 
published November 17, 1999 and incorporated by this reference, an estimate of tiie 
clipping error is generated at the receiver, and is used to reconstinict a frequency domain 
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compensation signal that is applied to the received signal, to remove the effects of any- 
clipping. 

[0013] A method of reducing the PAR in DMT transmissions without involving 

a loss of data rate is described in Gatherer and PoUey, "Controlling cHpping probabiUty 

5 in DMT transmission". Proceedings of the Asilomar Conjerence on Signals, Systems, and 
Computers, (1997), pp. 578-584, incorporated herein by this reference. As noted above, 
one function carried out in the training sequence on initiatian of a DSL session 
determines the number of bits per symbol (i.e., the bit loading) assigned to each 
subchannel. After bit loading, it is common for a number of subchannels (t5rpically at 

10 higher frequencies) to remain xmloaded, carrying no data symbols at all. In the Gatherer 
and Polley article, the PAR is reduced by using these unloaded subchannels to carry a 
"signal" that contains no payload, but that has the effect of reducing the ampUtude of 
the time domain signal to below the PAR amplitude threshold, in most cases. 

[0014] In summary, the Gatherer and Polley method performs an iterative 

15 process to derive the symbols for the unloaded subchannels. In short, an initial trial 
value (possibly zero) of the unloaded subchannel signal is added, in the time-domain, to 
the time-domain signal after the IDFT. A nonlinear function corresponding to the 
cUppkig amplifier is applied to the summed signal; if no dipping results (i.e., if none of 
the signal elements change), the current trial value of the unloaded subchannel signal is 
20 kept as part of the signal. If, on the other hand, the nonlinear clipping amplifier fimction 
indicates dipping, Ihe clipping is used to determine a new trial signal for the xmloaded 
subchannels, and the process repeated until clipping does not occur. 

[0015] Referring now to Figure la, an example of transmission of a DMT signal 

over a DSL communication system according to the Gatherer and Polley artide will now 
25 be described. Transmitting modem 10 receives an input bitstream that is to be 
transmitted to a receiving modem over a transmission channel This input bitstream is a 
serial stream of binary digits, in the appropriate format as produced by the data source, 
and is received by bit-to-symbol encoder 11 in transmitting modem 10. Encoder 11 
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groups the bits in the input bitstream into multiple-bit symbols that are used to 
modulate the DMT subchannels. The number of bits in each subchannel varies 
according to the bit loading assigned to each subchannel in the initialization of the 
conmiunication session, as known in the art. In addition, encoder 11 may also use error 
5 correction coding, such as Reed-Solomon coding, for error detection and correction 
purposes; other types of coding, such as treUis , turbo, or LDPC coding, may also be 
applied for additional signal-to-noise ratio improvement. The symbols generated by 
encoder 11 are typically complex symbols, including both amphtude and phase 
information, and correspond to points in the appropriate modttlation constellation (e.g., 
10 quadrature ampKtude modulation, or QAM). 

[0016] The encoded symbols are then applied to inverse Discrete Fourier 

Transform (IDFT) (also referred to as inverse Fast Fourier Transform (IFFT)), fxmction 
12. IDFT function 12 associates each input S5nnbol with one subchannel in the 
transmission frequency band, and generates a corresponding number of time domain 

15 symbol samples according to the Fourier transform. These time domain symbol samples 
are then converted into a serial stream of samples by parallel-to-serial converter 13. 
Functions 11 through 13 thus convert the input digital bitstream into a serial sequence of 
symbol values representative of the sum of a number of modulated subchannel carrier 
frequencies, the modulation indicative of the various data values. Typically, N/2 

20 unique complex symbols (and its N/2 conjugate symmetric symbols) in the frequency 
domain will be transformed by IDFT function 12 into a block of N real-valued time 
domain samples. 

[0017] As known in the art, function 14 adds a cyclic prefix to each block of 

serial samples presented by parallel-to-serial converter 13, by copying a selected number 
25 of sample values from the end of the block, and prepending the copy to the beginning of 
the block. The cyclic prefix has the effect of limiting intersymbol interference (ISI) due to 
energy from a previous symbol spreading into the next symbol due to the channel 
response. The cyclic prefix causes the datastream to appear to be periodic over a block 
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of N of the N+P samples, where P is the length of the prefix, so that the equivalence 
between frequency domain multiplication and time-domain convolution is valid. 

[00181 Upsampling function 15, and digital filter 16, then process the digital 

datastream in the conventional manner. As known in the art, upsampling function 
5 doubles or quadruples (or applies any multiple) tiie datasfaream to increase tiie sample 
rate, by inserting zero-value samples between each actiial signal sample. Digital filter 16 
may include such operations as a digital low pass fUter for removing image components, 
and digital high pass filtering to eliminate POTS-band or ISDN interference. The 
digitally-filtered datastream signal is then converted into the analog domain, by digital- 
la, 10 to-analog converter 17. Analog filtering (not shown) may then be performed on tiie 
5 output analog signal, such filtering typically including at least a low-pass filter. The 

W analog signal is then amplified by amplifier 18 which, according to this embodiment of 

asi 

% the invention, includes a clipping function, such as a hardlimiting cHpping function that 

W limits the ampUtiide to a maximum (both positive and negative polarities). As described 

r 15 in U.S. Patent No. 6,226^22, digital-to-analog converter 17, amplifier 18, and any analog 

S filtering, may be implemented in coder/ decoder (codec) integrated circuit in 

fiJ transmitting modem 10 (Figure 2) . 

S| 

9 [0019] The amplified analog output is then applied to a transmission channel, 

for forwarding to a receiving modem. According to conventional ADSL technology, the 
20 tiansmission channel consists of some length of conventional twisted-pair wires. In 
general, the receiving modem (not shown) reverses the processes performed by 
tiransmitidng modem 10 to recover the input bitstiream as tiie tiransmitted 
commtmication. 

[0020] An example of the hardlimiting cUpping function appUed by amplifier 18 

25 in this conventional transmission system is illustrated in Figure lb. Considering the 
input to amplifier 18 as sampled analog signal Xs(n), amplifier 18 generates an output 
signal corresponding to the function /(Xs(n)), as shown in Figure lb. This function 
/(Xs(n)) can effectively be expressed as: 



in 



is 



ytxs(n)) = 



-M, forBx,(n)<-M 
Bx,(n), for|Bx3(n)|<M 
M, forBx3(n)>-M 



where B is the gain of amplifier 18, and where M is the maxiinmn output magnitude of 
amplifier 18, as shown in Figure lb. As discussed above, each time tiie amplitude of a 
sample of signal Xs(n) exceeds the clipping threshold | Bxs(n) | , data may be lost from the 
5 payload signal; alternatively, if the average power is held to be low enough to effectively 
avoid dipping, amplifier 18 will not be operating efficiently, and resolution in the 
transmitted signal will also be compromised. 

[0021] According to tiie Gatherer and Polley metiiod, tmloaded channel 

encoding function 19 generates symbols that are assigned to unloaded subchannels in 
W 10 order to reduce the likelihood of clipping, by reducing the peak-to-average ratio (PAR) 

SS5! . 

of the transmitted signal. A cUp prevention signal Xc can be derived that is orthogonal 
to frequency-domain payload signal Xs, by assigning dip prevention signal Xc to 
subchannels where payload signal Xs is zero. An indicator matiix G can thus be derived 
as a diagonal matrix with ones corresponding to unloaded subchannels, and zeroes 
15 elsewhere. This indicator matrbc G will thus have the properties: 



GXs=0 
GXc ~ Xc 

indicating that signals Xc and Xs are orthogonal to one anotiier. 

[0022] To reduce the peak amplitude of payload Xs after IDFT and D/A 

20 conversion, cUp prevention signal Xc is selected to maintain payload signal Xs below the 
clipping threshold: 

ytxs + xc)-(xs + xc)=0 

in the time domain, where Xs = pHXs, pH being the inverse DFT operator (and F being the 
DFT operator). In the frequency domain, one may rewrite this relationship as: 
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ytFHXs + FHXc) - (FHXs + FHXc) = 0 
One may readily solve for dip prevention signal Xc in the frequenqr domain by: 

Xc = Fy(FHXs + FHXc)-Xs. 

The ortiiogonality constraint between payload signal Xs and clip prevention signal Xc is 
5 enforced using the indicator matrix G: 

Xc = GXc = G Fy^FHXs + FHXc) - GXs 
=GFy?FHXs*FHXc) 

According to the Gatherer and Polley article, the cHp prevention signal Xc is solved for 
using the Projection Onto Convex Sets (POCS) method. For example, an iterative method 
10 can be used to solve for the clip prevention signal Xc 

1. Start witii an initial trial value for dip prevention signal Xc (e.g., Xc=0); 

2. Form the time domain signal )fe + Xc = FHXs + FHXc; 

3. Apply the nonlinear dipping function to the time domain signal Xs + Xc 
and determine whether any elements change (i.e., are dipped). If not, 

15 the current faial value of dip prevention signal Xc is adequate. 

4. If dipping occurred, perform a DFT of tiie dipped time domain signal 
Xs + Xc and form a new b:ial value for clip prevention signal Xc using the 
above equation induding flie indicator function, and repeat the process. 

An alternative approach, according to the Gatherer and PoUey method, is to perform tiie 
20 step of forming a new trial value for the dip prevention signal Xs in the time domain, using 
a vector of preconqjuted values and no additional tiansforms. This is accomplished by 
performing an IDFT of the expression for the dip prevention signal Xc, and then 
precalculating a shaping matrix S=FhGF that is appHed to tiie dips of the dipped time 
domain signal Xs + Xc to be added to tiie prior tiial value of dip prevention signal Xc to 
25 produce tiie next trial value. The shaping mateix need only be calculated once for a given 
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indicator matrix G, and the time domain transform Xs=FhXs need only be performed once, 
according to this time-domain update approach. In this approach, the resulting time- 
domain dip prevention signal Xc is added in after IDFT function 12, as shown in Figure la. 

[0023] A partictilar benefit of the Gatherer and Polley method is that die receiver 

need not transform the signal, and that no data rate is lost. The modifying signal is 
applied by the transmitter only affects unloaded subchannels, and thus is not considered 
by the receiver in demodulating the transmitted sigiwl. Indeed, titie payload portion in 
the loaded subchannels are not modified in any way. 
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BRIEF SUMMARY OF THE INVENTION 

[0024] It is an object of this invention to provide a method and system for 

reducing the peak-to-average ratio (PAR) of a discrete multitone (DMT) modulated 
signal in an improved manner. 

5 10025] It is a further object of this invention to provide such a method and 

system in which ttie effects of upsampling and filtering are accounted for in reducing the 
PAR. 

[0026] It is a further object of this invention to provide such a method in which 

the PAR reduction can be efficiently accomplished by the execution of a software routine 
10 by programmable logic, such as a digital signal processor (DSP). 

[0027] Other objects and advantages of the present invention will be apparent to 

those of ordinary skiU in the art having reference to tiie following specification together 
with its drawings. 



[0028] The present invention may be implemented by an iterative approach 

3 15 executed in the transmitter loop of a DMT communications system, such as a DSL 
modem, in which upsampling and filtering is applied after Inverse Discrete Fourier 
Transform (IDFT) modulation of the payload symbols. Subchaimels are identified that, 
in the training process, were loaded with no symbols. An initial trial value of an 
orthogonal cUp prevention signal is selected; the initial value may be zero, and the clip 
prevention signal is orthogonal in the sense that it is appUed to the unloaded 
subchannels. A dip flag is reset, and a time domain sum signal is produced by 
performing the IDFT of the payload signal and the trial value of the dip prevention 
signal, for the first phase of the polyphase filter corresponding to the upsampling and 
filter, and summing the results. A nonlinearity corresponding to the cUpping threshold 
25 is appHed to the summed signal, and the result is examined to determine if any cUpping 
resulted; if so, the cHp flag is set. A new estimate for the dip prevention signal is then 
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generated, and the IDFT and nonlinear clipping function is applied for the next phase of 
the polyphase filter. Upon completing the process for each phase of the polyphase filter, 
the clip flag is interrogated. If set, the process is repeated; if not, the most recent trial 
value of the clip prevention signal is maintained as its most recent trial value, and is 
5 transmitted in the imloaded channels along with the DMT modulated payload. 
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BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING 

[0029] Figure la is a functional block diagram of the functions involved in 

transmission of a DMT modulated signal, according to the prior art. 

[0030] Figure lb is a transfer characteristic of a conventional clipping amplifier. 

5 [0031] Figure 2 is an electrical diagram, in block form, of a transceiver 

constructed according to the preferred embodiments of the invention. 

[0032] Figure 3 is a functional block diagram of communications between two 

DMT transceivers, where at least the transmitting transceiver is constructed according to 
a first preferred embodiment of the invention. 

10 [0033] Figure 4 is a polyphase representation of ihe operation of the upsampling 

and filtering operations performed in the transmitting transceiver according to the 
preferred embodiments of the invention. 

[0034] Figure 5 is a flow diagram illustrating the operation of a first preferred 

embodiment of the invention in generating a clip prevention signal. 

15 [0035] Figure 6 is a functional block diagram of communications between two 

DMT transceivers, where at least the transmitting transceiver is constructed according to 
a second preferred embodiment of the invention. 

[0036] Figure 7 is a flow diagram illustrating the operation of the second 

preferred embodiment of the invention in generating a cUp prevention signal. 

20 [0037] Figure 8 is a flow diagram illustrating the process of updating a trial 

value of the dip prevention signal according to the second preferred embodiment of the 
invention as shown in Figure 7. 
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DETAILED DESCRIFHON OF THE INVENTION 
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[0038] The present invention will be described in connection with an example of 

its implementation in a commtmications transceiver, such as a Digital Subscriber Line 
(DSL) modem. It will be apparent to those skilled in the art having reference to this 
5 specification that this invention is particularly well-suited for use in such an application. 
However, it is also contemplated that this invention will be of similar benefit in many 
other applications that involve multicarrier modulation. It is therefore to be understood 
that these and other alternatives to the embodiment described below are contemplated 
to be within the scope of the invention as claimed. 

10 [0039] Referring now to Figure 2, the construction of DSL modem 100, in the 

form of a transceiver (i.e., for both transmitting and receiving) and constructed 
according to the preferred embodiments of the invention, will be described. The 
particular architecture illustrated in Figure 2 is presented by way of example, as those 
skilled in the art having reference to this specification will readily appreciate that this 

15 invention may be implemented into DSL modems of other architectures, and also into 
other communications devices for transmitting DMT modulated signals. In particular, 
while particular functions will be described as being performed by certain integrated 
circuit devices in this description of modem 100, it will be appreciated by those skilled in 
the art that alternative arrangements of the particular integrated circuits within modem 

20 100 may be used in connection with this invention, for example by integrating more 
functions to within a fewer number of integrated circuits. In addition, modem 100 may 
correspond to one port of a central office modem, which may of course carry out 
simultaneous transmissions with many client modems, or alternatively modem 100 may 
correspond to a client modem for commtmicating with a single central office for DSL 

25 access. 

[0040] As shown in Figure 2, modem 100 includes hybrid circuit 30, which is 

coimected to twisted-pair transmission channel H. Hybrid circuit 30 is a conventional 
circuit that converts the two-wire arrangement of twisted-pair f aciUty H to dedicated 
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transmit and receive lines connected to line driver and receiver 32, considering that 
modem 100 is suitable for carrying out both the receipt and the transmission of DMT 
modulated communications. Line driver and receiver 32 is a high-speed line driver and 
receiver for driving and receiving ADSL signals over twdsted-pair faciUty H; an example 

5 of a suitable line driver and receiver circuit is the THS7102 line driver/receiver available 
from Texas Instruments Incorporated. Line driver and receiver 32 is bidirectionally 
coupled to coder/ decoder (codec) circuit 34, which carries out analog and digital 
filtering, and analog-to-digital conversion. An example of a suitable codec device may 
be selected from the TLV320ADlx device family available from Texas Instruments 

10 Incorporated. Hybrid circuit 30, line driver and receiver 32, and codec 34 are often 
referred to, in the aggregate, as an "analog front end". 

[0041] Modem 100 also includes DSP 35, which serves as a digital transceiver for 

DSL communications. As noted above, DSP 35 is a higji performance digital signal 
processor, for carrying out the digital operations of the present invention as will be 

15 described below; exemplary devices suitable for use as DSP 35 include the TMS320C6x 
family of DSPs available from Texas Instruments Incorporated. For the example in 
which modem 100 is a client modem, conventional interface 36 couples DSP 35 in the 
usual manner to a data source, which in this example is cUent workstation W; 
alternatively, if modem 100 is implemented at a central office, the data source would be 

20 the appropriate communications network or switch at the central office. 

[0042] Figure 3 ftmctionaUy illustrates an example of a DSL communication 

system into which the preferred embodiment of the invention is implemented. The 
same reference numerals will be used in Figure 3 for the same functions performed in 
conventional modem 10 as described above relative to Figure 1. In the system of Figure 
25 3, only one direction of transmission (from transmitting modem 100 to receiving modem 
200) is iUustrated; it will of course be understood by those skilled in the art that data will 
also be communicated in the opposite direction (from modem 200 to modem 100). 
According to this embodiment of the invention, the DSL communications carried out in 
the system of Figure 3 are of the asymmetric DSL (ADSL) type, la this arrangement, so- 
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called "downstream" communication is from the telephone company central office to 
the subscriber's modem while "upstream" communication is from the subscriber to the 
central office. 

[0043] As shown in Figure 3, transmitting modem 100 receives an input 

bitstream that is to be transmitted to receiving modem 200. The input bitstream may be 
gaierated by a computer at the same location as transmitting modem 100, or 
alternatively and more likely for the case of a central office modem, is communicated 
over a computer network, in the Internet sense, that is coupled to transmitting modem 
100. The input bitstream is a serial stream of binary digits, in the appropriate format as 
produced by the data source. 

[0044] The input bitstream is received by bit-to-symbol encoder 11 in 

transmitting modem 100. Encoder 11 performs multiple functions in this embodiment 
One such function is to group the bits in the input bitstream into multiple-bit symbols, 
for example ranging up to as many as fifteen bits. The symbols will be the values used 
to modulate the various subchannels, as wiU be described below. The symbols typically 
are not necessarily representative of the transmitted data itself, but instead simply 
correspond to grouping of the serial bitstream for broadband transmission. The number 
of bits in each symbol produced by encoder 11 wiU vary, in this DSL appUcation, 
according to the bit loading assigned to each subchannel in the initialization of the 
communication session, as known in the art. hi addition, encoder 11 may also use error 
correction coding, such as Reed-Solomon coding, for error detection and correction 
purposes; other types of coding, such as trelKs, turbo, or LDPC coding, may also be 
appUed for additional signal-to-noise ratio improvement. In addition, the symbols 
generated by encoder 11 are typicaUy complex symbols, including both amphtude and 
phase information, and correspond to points in the appropriate modulation constellation 
(e.g., quadrature amplitude modulation, or QAM). 

[0045] The encoded symbols are then appUed to inverse Discrete Fourier 

Transform (IDFT) (also referred to as inverse Fast Fourier Transform (IFFT)) function 12. 
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IDFT fxinction 12 associates each input symbol with one subchannel in the transmission 
frequency band, and generates a corresponding number of time domain symbol samples 
according to the Fourier transform. These time domain symbol samples are then 
converted into a serial stream of samples by parallel-to-serial converter 13. Functions 11 
5 liirough 13 thus convert the input digital bitstream into a serial sequence of symbol 
values representative of the sum of a number of modulated subchaimel carrier 
frequencies, the modulation indicative of the various data values. Typically, N/2 
imique complex symbols (and its N/2 conjugate synunetric symbols) in the frequency 
domain wiU be transformed by IDFT function 12 into a block of N real-valued time 
10 domain samples. 

[0046] In function 14, a cyclic prefix is then added to each block of samples 

presented by IDFT function 12 through parallel-to-serial converter 13. The cyclic prefix 
is constructed by a selected number of sample values at the end of the block, and these 
values are then prepended to the block by function 14. For example, if each block 

15 includes 512 samples, a typical cyclic prefix is on the order of thirty-two samples. The 
cyclic prefix has the effect of limiting inters5nnbol interference (BI) due to energy from a 
previous symbol spreading into the next symbol due to the channel response; in other 
words, the cycUc prefix causes the datastream to appear to be periodic, so that the 
equivalence between frequency domain multiplication and time-domain convolution is 

20 valid over N of the N+F samples (F being the number of samples in the prefix). Those 
skilled in the art having reference to this specification wiU readily recognize that each of 
functions 11 through 14 may be carried out, and preferably are carried out, as digital 
operations. As such, it is contemplated tiiat a high-performance digital integrated 
circuit, such as a digital signal processor (DSP) having computational power similar to 

25 or greater than the TMS320C5x and TMS320C6x DSPs available from Texas Instruments 
Incorporated will be suitable for these functions. 

[0047] Upsampling function 15, and digital filter 16, then process the digital 

datastream in the conventional manner. As known in the art, upsampling function 
doubles or quadruples (or applies any multiple) the datastream, by inserting zero- 



valued samples between each signal sample, thus increasing the sample rate. Digital 
filter 16 may include such operations as a digital low pass filter for removing image 
components, and digital high pass filtering to eliminate POTS-band or ISDN 
interference. The digitally-filtered datastream signal is then converted into the analog 
domain, by digital-to-analog converter 17. Analog filtering (not shown) may then be 
performed on the output analog signal, such filtering typically including at least a low- 
pass filter. The analog signal is then amplified by amplifier 18 which, according to this 
embodiment of the invention, includes a clipping fimction, such as a hardlimiting 
clipping fxmction that limits the amplitude to a maximum (both positive and negative 
polarities). As described in U.S. Patent No. 6,226,322, digital-to-analog converter 17, 
amplifier 18, and any analog filtering, may be implemented in coder/ decoder (codec) 
integrated circuit in transmitting modem 100 (Figure 2). 

[0048] The amplified analog output is then applied, via the appropriate line 

driver and other drcuitry (not shown in Figure 2) to transmission channel H, for 
forwarding to receiving modem 200. According to conventional ADSL technology, the 
transmission chaimel consists of some length of conventional twisted-pair wires. The 
transmission channel H communicates the transmitted analog signal to receiving 
modem 200, which, in general, reverses the processes performed by transmitting modem 
100 to recover the input bitstream as the output bitstream. 

[0049] As shown in Figure 3, the functions of modem 200 in processing received 

downstream communications according to the preferred embodiment of the invention 
will now be described, according to an example based upon time domain equalization. 
It is of course tmderstood that the recovery of the payload signal at the receiver may be 
alternatively carried out according to other techniques; because the present invention 
does not require any change at the receiver, such alternative receiver approaches may 
also be used in connection with this invention. In the example of Figure 3, following 
physical receipt of the signal through hybrid circuit 30 and line driver and receiver 32 
(Figure 2), analog filtering function 21 is performed by codec 34, primarily to remove 
high frequency noise and aliasing from the incoming signal; as described in U.S. Patent 
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No. 6,226,322, commonly assigned herewith and incorporated by this reference, 
equalization of the analog signal may also be performed to compensate for line 
attenuation characteristics of transmission channel H. Analog-to-digital conversion 
process 22 then converts the filtered analog signal into the digital domain, foUowing 
which conventional digital filtering function 23 is appUed to augment the function of the 
analog filters. Functions 21, 22, and 23 are preferably performed by codec device 34 of 
Figure 2, to provide a filtered digital data stream to DSP 35 for digital processing, as will 
now be described. 

[0050] Digital filter function 23 forwards the filtered digital datastream to time 

domain equalizer (TEQ) 24. TEQ 24 is preferably a finite impulse response (FIR) digital 
filter, implemented as a software routine executed by DSP 35 (of Figure 2), that is 
designed to effectively shorten the length of the impulse response of the transmission 
channel H, including the filtering that is performed prior to receipt by TEQ 24. The 
design of this TEQ filter is realized by the selection of the particular coefficients of the 
FIR implementing TEQ function 24 during initialization, or "training" of modem 200 in 
combination with modem 100 upon establishing a communications session. 

[0051] Following TEQ 24, the cycEc prefix is removed from each received block 

in function 25. Serial-to-parallel converter 26 then converts the datastream into a 
number of samples (2N) for application to Discrete Fourier Transform (DFT) function 27. 
Because the received datastream appUed to DFT function 27 is a time-domain 
superposition of the modulated subchannels, the DFT of this datastream wiU recover the 
modulating symbols at each of the subchannel frequencies, reversing the IDFT 
performed by function 12 in the transmission processing by transmitting modem 100. 
The output of DFT function 27 is a frequency domain representation of the transmitted 
symbols multipUed by the frequency-domain response of the effective transmission 
channel, under the assumption that liie cyclic prefix is at least as long as the (possibly 
shortened) channel response h(t). Frequency-domain equalization (FEQ) function 28 
then divides out the frequency-domain response of the effective channel, recovering the 
modulating symbols. Symbol-to-bit decoder function 29 then resequences the symbols 
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into a serial bitstream, decoding any encoding that was applied in the transmission of 
the signal, and producing an output bitstream that corresponds to the input bitstream 
upon which the transmission was based. This output bitstream is then forwarded to the 
client workstation, or to the central office network, depending on the location of 
5 receiving modem 200. 

[00521 According to this embodiment of the invention, transmitting modem 100 

also includes unloaded subchannel encoding function 119. According to the preferred 
embodiments of the invention, encoding function 119 effectively inserts symbols into 
one or more subchannels prior to IDFT 12, corresponding to one or more subchannels 
10 that have a bit loading of zero determined in the training routine. As such, the operation 
of function 119 adds information into the subchannels that are not otherwise loaded 
with payload data from the input bitstream. According to this preferred embodiment of 
the invention, however, this information is generated in such a manner as to incorporate 
the effects of upsampling function 15 and filter 16, as will now be described. 

15 [0053] As known in the art, upsampling and filtering as performed by 

upsampling function 15 and filter 16 in order to increase the sample rate, and eliminate 
image components and POTS-band or ISDN interference. The unloaded channel 
encoding described in the Gatherer and Polley article, while generaUy usefiil, does not 
take this upsampling and filtering into account. As a result, time-domain samples that 

20 do not exceed the cHpping threshold prior to upsampling and filtering, and therefore for 
which a compensating cHp prevention signal Xc would not be generated, may exceed the 
dipping threshold after upsampling and filtering but without having encoded unloaded 
subchannels. Because of the effects of upsampling and filtering, the effects of encoding 
the unloaded subchannels as described in the Gatherer and Polley article are limited in 

25 some applications. 

[0054] The theory of operaticm of ihe present invention will now be described, 

first in connection with filtering in the absence of upsampling, and then in connection 
with both upsampling and filtering prior to digital-to-analog conversion. The goal of 
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the method according to this invention is to derive a clip prevention signal Xc (in the 
frequency domain; Xc in the time domain) such that its sum with the payload signal Xs 
does not clip after digital filter fimction 16 is applied. 

10055] Consider the impulse response of filter fxmction 16 as /i (and upsampling 

function 15 to upsample by a factor of one, in other words adding no samples). In the 
time domain, the goal of encoding function 119 is to apply a frequency domain dip 
prevention signal Xc that has a time domain transform Xc for which: 

y(/i*Xs + h*Xc) - (fe*Xs + /l*Xc) = 0 
where * is the time-domain convolution operator. In the frequency domain: 

y(FHHXs + FHHXc) - (FHHXs + FHHXc) =0 

where H is a diagonal matrix having non-zero entries corresponding to the DFT of the 
impulse response ft. To the extent that interblock effects from filter function 16 are 
important, these effects could be handled by saving a previous state of the filter, or by way 
of a conventional overlap and add or save method. To define the dip prevention signal Xc 
one may rearrange this frequency domain equation: 

FHHXc =y(FHHXs + FHHXc) " pHRXs 
and appl5nng the DFT operator F 

HXc = F^FHHXs + FHHXc) - HXs 

To isolate dip prevention signal Xc one may multiple by the inverse H-i of the filter 
response H: 

Xc = H-iFytpHHXs + FHHXc) - Xs 

Recalling the orthogonality constraint, using indicator matrix G which is a diagonal matrix 
with ones at the entries having unloaded subchannels: 
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Xc GXc 

= GH-^F/FHHXs + FHHXc) - GXs 
and because GXs = 0 from the orthogonality: 

Xc = GH-iFy^FHHXs + F^HXc) 

5 This fixed point equation for clip prevention signal Xc can now be solved by way of 
Projection Onto Convex Sets (POCS). One should note that, because clip prevention signal 
Xc is multipUed by the filter response H, if this filter response H is small on unloaded 
subchannels, then the resulting clip prevention signal Xc wiU likely be ineffective. 

[0056] Referring to Figure 3, if one considers the time domain input to 

W 10 upsampling function 15 as signal x(n), and the output of filter function 16 as y(n), with 

an upsampling factor of L, an equivalent polyphase representation is illustrated in 
f^' Figure 4, As shown in Figure 4, input signal x(n) is effectively applied to L filter 

? fimctions 16o through 16l-i, with corresponding filter responses ho through hi-i, with the 

n I output of each filter function 16j represented as signal yj(n). The individual polyphase 

^ 15 filters h for 1 from 0 to L-1 are related to the original filter response h by: 

P fc(n)=ft(Ln+0 

The filtered signals yo(n) through yL-i(n) are then each upsampled by a factor L by 
upsampling functions 15. The upsampled outputs are combined into output signal y(n) by 
adder 40, with the higher order upsampled signals yi(n) through yL-i(n) delayed by 1 
20 through L-1 delay stages 42. If the individxial filter output signals y;(n) do not clip, the 
combined output y(n) will not dip. This is because upsampling functions 15 each insert 
zero-valued samples between each sample of individual output signals yz(n), and the delay 
stages 42 delay the mxiltiple paths so that the filtered output signals yi(n) do not add to one 
another. 

25 [0057] Recalling the fixed-point relationship of dip prevention signal Xc: 
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Xc = GH-iFy(FHHXs + FHHXc) 
the polyphase analysis indicates that the upsampling factor L may be incorporated by: 

Xc = GHfi FytTOXs + FhH/Xc) 

for Z = 0, 1, . , . , L-1, where Hi is the frequency response of the polyphase filter 16/. This 
solution for clip prevention signal Xc thus satisfies the L polyphase equations. 

[0058] Based on this theory of operation and referring now to Figure 5, a method 

of performing encoding fxmction 119, by determining the clip prevention signal Xc, 
according to the preferred embodiment of the invention wiU now be described in detail. 
As noted above, encoding function 119 is contemplated to be implemented by way of 
DSP 36 (Figure 2), with the method steps shown in Figure 5 carried out by DSP 36 
executing a software routine. Of course, other realizations of function 119 are also 
contemplated to be within the scope of this invention, such realizations including 
custom logic, for example in the form of a coprocessor, residing within modem 100. 

[0059] As shown in Figure 5, the method begins with process 50, in which the bit 

loading of ihe DMT subchannels is determined, and those subchannels receiving no bit 
loading are identified. It is contemplated that process 50 will be carried out in the 
initialization of the DSL session, by way of a conventional modem "training^' routine; it 
is also contemplated that the bit loading may be updated periodically during the DSL 
session, to account for changes in the transmission environment. In any case, certain 
subchannels within the transmission bandwidth are identified in process 50 as being 
unloaded, and an indicator matrix G is defined. As noted above, indicator matrix G is a 
diagonal matrix with ones at the locations corresponding to the xmloaded subchannels, 

[0060] In process 52, DSP 36 encodes a block of symbols from the input 

bitstream, which may be generated by a computer at the same location as transmitting 
modem 100, or received over a computer network. As discussed above relative to 
encoder function 11 of Figure 3, the input bitstream is a serial stream of binary digits 
that is grouped, in process 52, into multiple-bit s5rmbols, with the number of bits in each 
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symbol depending on the bit loading determined in process 50. Other processes that 
may optionally be carried out in connection with process 52 include error detection and 
correction coding, or trellis , turbo, or LDPC coding. The output of process 52 is 
contemplated to be sequences of complex symbols, including both amplitude and phase 
5 information, and corresponding to points in the modulatian constellation for the 
particular communications class (e.g., quadrature amplitude modulation, or QAM). In 
process 54, these symbols are associated with subchannels in the transmission band (i.e., 
those subchannels having non-zero bit loading), and together comprise payload signal 
X^(n). 

10 [00611 As noted above, the goal of the method of Figure 5 is to derive clip 

prevention signal Xc that, when summed with payload signal Xs and converted into the 
time domain, has amplitude peaks that are within the clipping threshold, including the 
effects of upsampling and fUtering. Because of the iterative nature of this method, 
process 56 begins by initializing a trial value for cHp prevention signal Xd it is 

15 contemplated that an initial tiial value may be, and perhaps preferably is, zero at all 
subchannels (including the loaded subchannels). Flag CLIP is reset in process 58, and a 
loop counter I is initialized to zero. Loop counter I counts the number of passes through 
the remainder of the method of Figure 6, and corresponds to a current phase of the 
polyphase representation of the L phases in Figure 4. 

20 [0062] In process 62, DSP 36 forms a time domain signal from the frequency- 

domain signal block Xs and the current trial value of clip prevention signal Xc. As 
discussed above, this time-domain signal is based upon a IDFT (or IFFT, as the case may 
be) of signals Xs, Xc, contemplating the effects of tiie current phase of the filter. This time 
domain signal thxis is generated by summing the results of the IDFTs: 

25 FHHjXs + FHHiXc 

Of course, for the trivial case of initial trial cHp prevention signal Xc=0, the second term 
of this sum is zero. As before, the Hj term is the frequency response of filter function 16i 
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for the current phase / imder investigation. The restilting time domain signal will, of 
course, be represented by a sequence of sample amplitudes. 

[0063] In process 64, the amplifier function / is applied to this time domain 

signal by DSP 36. While it is contemplated that this function /is typically a hardUmiting 
5 clipping amplifier as shown in Figure lb, the particular function /applied in process 64 
may correspond to any linear or nonlinear function (but typically including 
nonlinearities). Decision 65 is then performed to determine whether any of the samples 
in the current time domain sum signal is cUpped by the application of function / in. 
process 64. If not (decision 65 is NO), decision 69 is tihen executed to determine whether 
10 loop cotmter I indicates that the fiaal phase has been evaluated. If not (decision 69 is 
NO), loop coxmter I is incremented in process 70, and a new time domain sum signal is 
produced in process 62 and the process repeated for the next phase. 

[0064] If clipping is detected (decision 65 is YES), flag CLIP is set in process 66 to 

indicate that a clipped sample has been detected. In process 68, a new trial value of clip 

15 prevention signal Xc is generated, by applying the operator GHr^ to the clipped time 
domain signal (output by process 64). This operator corresponds to the inverse of the 
filter response for the current phase, multiplied by the indicator matrix G, and produces 
a new trial value for clip prevention signal Xc, which will be used in the next iterations 
through the process. Decision 69 is then repeated to determine whether additional 

20 phases remain to be analyzed. If so (decision 69 is YES), loop counter / is incremented in 
process 70, and the process repeated for the next phase using the new trial clip 
prevention signal Xc beginning from process 62. 

[0065] Upon decision 69 indicating (with a YES result) that loop coimter I has 

reached its terminal value of L-1 (the value L being the upsampling ratio), decision 71 
25 interrogates flag CLIP. If flag CLIP is set (decision 71 is YES), one of the time-domain 
signals for one of the phases was clipped, and a new trial value for clip prevention signal 
Xc is to be generated, for at least one phase within the most recent loop through 
processes 62 et seq. In this case, control passes to process 58, in which flag CLIP is reset, 
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and to process 60, in which loop counter I is reinitialized to zero, and the process is 
repeated for the same input signal Xs but with the new trial dip prevention signal Xc, 
beginning with the zeroeth phase. 

[0066] Once the process has been repeated for all L phases (decision 69 is YES) 

and witii none of the phases having been clipped by application of the function /ifci 
process 64 (decision 71 is NO), Ihe current trial value of clip prevention signal Xc can be 
deemed to be effective in eliminating dipping from the payload signal Xs. As described 
above, this generation of clip prevention signal Xc is performed in a manner that takes 
into accoimt upsampling fimction 15 and filter 16. Process 72 is then performed, by way 
of which DSP 36 loads the current imloaded subchannels with the current clip 
prevention signal Xc. This new signal is then passed to IDFT function 12, and the 
modulation and transmission of the corresponding modulated signal is carried out. 

[0067] According to this embodiment of the invention, therefore, the peak to 

average ratio (PAR) is reduced by the generation of a dip prevention signal, and by the 
encoding of otherwise tmloaded subchannels with that dip prevention signal. Because 
the dip prevention signal is orthogonal to the payload signal, no loss of data is involved 
in this PAR reduction. In addition, because the clip prevention signal resides 
exclusively in unloaded subchannels, the clip prevention signal is completely 
transparent to receiving modem (e.g., modem 200 ia Figure 3) and involves no 
processing on the receiver end to determine whether clipping occurred or to compensate 
for any clipping. Rather, the subchannels carrying the cUp prevention signal are simply 
ignored, because they are unloaded subchannels insofar as the receiver is concerned. 
This method is also quite computationally efficient, considering that the only somewhat 
complex step involves the IDprs used in generating the time domain signal (process 
62). 

[0068] An altemative approach for deriving the clip prevention signal will now 

be described. As wiU become apparent from the following description, this approach 
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operates as a time domain update, and is contemplated to be a higher performance 
method, by using a vector of precomputed values and no additional transforms. 

[0069] Recalling the expression for clip prevention signal Xci 

Xc = GHr^Fy(FHH/Xs + FhHzXc) 
5 One can multiply both sides of this equation by the filter response Hi for the current phase: 

H/Xc = GF/(FHH/Xs + pHHzXc) 

= GFy^fe*Xs + /iz*Xc) 
^ The IDFT of both sides of this equation iato the time domain provides: 

y fti*Xc = FHGF/t/ii*Xs + /z/*Xc) 

€1 10 or, defining a shaping function S = F^^GF: 

f flf * Xc = S/lhi * Xs + fei * Xc) 

III This can now be rewritten to apply the update of the clip prevention signal by way of 

applying the shaping function to the clip updating: 



ft/ * Xc = S/{hi * Xs + fe/ * Xc) + S(ft/ * Xs + A/ * Xc) - S(/t/ * Xs + ft; * Xc) 
15 =S(/f;fe*Xs + /l/*Xc)-(fe*Xs + ftl*Xc))+S(fc*Xs + /?/*Xc) 

= S(/{hl * Xs + * Xc) + (ft/ * Xs + * Xc)) + fc * Xc 

Any perturbation in the 1^ must be translated to the other phases, however. Referring back 
to the frequency domain equation: 

HXc = GFytft/ * Xs + fti * Xc) 
20 this can be similarly rewritten as: 

HiXc = GFytfti*Xs + fti*Xc)+GF(ft?*Xs + ft/*Xc)-GF(ftz*Xs + ft/*Xc) 
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= GF(4fc * Xs + ftz * xc) - (fe * Xs + fci * Xc)) + GFQii * Xs + ftr * Xc) 

= GF (/(/z/ * Xs + /Z/ Hc Xc) - (/Zi * Xs + hi * Xc)) + HiXc 
By midiiplying each side of this equation by HmH"^: 

HznXc = HxnHfiGF (/(/Zi * Xs + * Xc) - (fe * Xs + * Xc) ) + HmXc 

5 The IDFT of this equation provides: 

h« * Xc = FHHmHriGF(Xft/ * Xs + A/ * Xc) - (/zr * Xs + /z/ * Xc)) + fem * 
or, for shaping function Sjm = : F^HmHriGF: 

p hm*Xc = Slm (Ah * Xs + /Zi * Xc) - (fc * Xs + /Z/ * Xc)) + /Zm * Xc 

y Shaping function S/m thus shapes the application of the clip from the phase to the update 

10 for the phase. 

[0070] Referring now to Figures 6 through 8, a method for effecting the encoding 

of unloaded subchannels in order to reduce the PAR according to a second preferred 
embodiment of the invention and following this time-domain update approach, will 
now be described in detail. As will be evident from the followiag description, this 
15 method is performed in a manner that takes upsampling function 15 and filter 16 into 
accoimt. 

[0071] Figure 6 illustrates the commtmications process in which the data is 

transmitted and received according to this second embodiment of the invention, with 
the same reference numerals being used for the same functions as described above 
20 relative to Figure 3. According to this second embodiment of the invention, however, 
encoding ftmction 219 produces time-domain clip prevention signal Xc, which is applied 
into the signal stream at a point in time after the IDFT of payload signal Xs is performed 
in IDFT function 12. Because of the orthogonality constraint of cHp prevention signal Xc 
relative to payload signal Xs, the subchannel frequencies corresponding to clip 
25 prevention signal Xc wiU be different from those of payload signal Xs, and wiU 
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correspond to the subchannel frequencies of the tinloaded channels. As shown in Figure 
6, for example, clip prevention signal Xc is applied to parallel-to-serial converter function 
13, by way of which time-domain clip prevention signal Xc is merged with time-domain 
payload signal Xs. 

5 [0072] The method according to this second embodiment of the invention 

begins, as before, with the determination of the bit loading of the DMT subchannels, and 
the identifying of those subchannels receiving no bit loading, in process 80. As before, 
process 80 is preferably performed during the initialization of the DSL session. Process 
82 then encodes a block of symbols from the input payload bitstream into multiple-bit 
= , 10 symbols according to the corresponding symbol constellation (e.g., QAM) and associates 
O these symbols with subchannels in the transmission band according to the bit loading of 

Ijj process 80, to produce payload signal Xs(n). 

[0073] In process 85, an IDFT routine is performed to transform the frequency 

domain payload signal Xs into time domain signal Xs. This conversion corresponds to 
15 IDFT function 12, except of course that the xmloaded subchannels have zero-valued 
components. In process 86, clip prevention signal Xc is initialized to a first trial value, 
preferably zero (i.e., Xc(n)=0 for aU n). Flag CLIP is reset in process 88, and loop coimter 
I is initialized to zero in process 90. 

[0074] In process 92, a time-domain sum is formed from the payload time 

20 domain signal Xs and the current trial clip prevention signal Xc, for the phase of the 
polyphase representation of the upsampling and filter functions 15, 16 as shown in 
Figure 4. This time domain signal is the simple sum of the signals Xs, Xc incorporating 
their time domain filters for that phase: 

* Xs + /iZ * Xc 

25 This time domain sum signal is then processed, in process 94, by ftmction / which for 
example may be a hardlimiting amplifier characteristic as discussed above. 
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[0075] Decision 95 is then performed to determine whether any samples in the 

stun time-domain signal had their ampUtude modified by cUpping process 94. If not, 
control passes to decision 99, in which the loop counter I is compared against the 
terminal value of L-1 to determine whether all L phases of the representation have been 
considered. If not (decision 99 is NO), loop counter I is incremented in process 100, and 
the process is repeated for tiie next phase, from process 92. 

[0076] If one or more sample amplitudes did change (decision 95 is YES), then 

some sample amplitiides of the sum of the payload signal Xs and the current ti:ial cUp 
prevention signal Xc exceeded the clipping threshold of the nonlinear function appHed in 
process 94, requiring an adjustment in cUp prevention signal Xc. Flag CLIP is set in 
process 96, indicating tiiat the dipping occurred. The farial value of cUp prevention 
signal Xc is then adjusted in process 98. As discussed above in connection with the theory 
of operation of this second preferred embodiment of the invention, in this time domain 
update approach, an adjustinent to cUp prevention signal Xs due to one phase of the 
polyphase representation wiH also affect the other phases. Adjustment process 98 is 
performed so tiiat these effects are also appUed to the other phases, as wiU now be 
described relative to Figure 8. 

[0077] In process 110, a vector R is determined for the current phase. Vector R is 

used repeatedly, for each phase in the polyphase representation, and tiierefore it is 
beneficial to calculate it only once, as performed in process 110 in ti:iis embodiment of 
the invention. According to this embodiment of the invention, vector R is derived as: 

R=/(fe*Xs + hl*Xc) - (/li*Xs + /li*Xc) 

Phase cotmt m is then initialized to zero in process 112. 

[0078] In process 114, the rafi' phase of clip prevention signal Xc is modified 

according to the cHpping effects from the Z* phase tiiat were detected in process 94 and 
decision 95. Recalling the shaping function Sim. 
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S/« = :FHHxnHriGF 

using indicator matrix G as determined in process 80, the adjustment to dip prevention 
signal Xc is updated for the phase, based on the clipping in the phase, by: 

which is performed in process 114. Decision 115 determines if clip prevention signal Xc is to 
be adjusted in any more phases, by interrogating whether count m has yet reached the 
terminal value of L-1. If not (decision 115 is NO), coimt m is incremented in process 116, 
and cfip prevention signal Xc is adjusted for the next phase. If so (decision 115 is YES), fhea 
clip prevention signal Xc has been adjusted in each of the phases, for the clipping effects 
upon the currait phase (Z), and control passes to decision 99, 

[0079] Decision 99 determines whether the loop count I has reached the terminal 

value of L-1, and thus determined whether the current trial value of dip prevention 
signal Xc has been evaluated for all phases. If so (dedsion 99 is YES), dedsion 101 then 
interrogates flag CLIP to determine if the trial value of clip prevention signal Xc caused 
any clipping. If so (decision 101 is YES), the current value of clip prevention signal Xc 
must now be evaluated, beginning with the resetting of flag CLIP, initializing of loop 
count Z, and repeating of the sequence of operations beginning with process 92. 

[0080] Upon an evaluation of a trial value of clip prevention signal Xc that does 

not result in any clipping (decision 101 is NO), the routine is complete. The current clip 
prevention signal Xc is thus adequate to avoid clipping by reduction of the PAR using 
unloaded subchannels, in which case the current sum time domain signal Xs + Xc can be 
forwarded to paraflel-to-serial conversion function 13 (Figure 3), from which the 
transmission process continues. 

[0081] This second embodiment of the invention, as in the first embodiment, 

provides the important advantages of reducing the PAR and thus reducing clipping of 
the transmitted signal, without any loss of data and while maintaining a reasonable 
average power so that downstream amplifiers can operate efficiently. It is also 
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contemplated that the time domain update approach of this second embodiment of the 
invention will provide additional efficiency and performance in many applications, 
considering that only a single IDFT operation need be performed for a given block of 
data, because the derivation and updating of the clip prevention signal is performed in 
5 the time domain, after the IDFT of the payload signal, with no repeated IDFT processes 
being performed. 

[0082] In addition, according to both of these disclosed embodiments of the 

invention, no operations are required on the part of the receiving modem, because orUy 
the unloaded subchannels are affected by the clip prevention signal, and these 
10 subchannels are necessarily ignored by the receiving modem. Accordingly, the PAR 

Q reduction accomplished by this invention can be performed at the transmitter, 

O 

y transparently to the receiving modem circuitry. 

J'j: [0083] The preferred embodiments of the invention, as described above, are 

directed to the modulation of signals according to the Fourier transform (i.e., the 

L. 15 IDFT/DFT, or IFFT/FFT, transform pairs), which is an invertible orthogonalizing 

JJ; transform. It is contemplated that this invention may also be applied to, and be of 

benefit in connection with, any orthogonalizing transform, many types of which are weU 

ff known in the art. 

[0084] While the present invention has been described according to its preferred 

20 embodiments, it is of course contemplated that modifications of, and alternatives to, 
these embodiments, such modifications and alternatives obtaining the advantages and 
benefits of this invention, wiU be apparent to those of ordinary skill in the art having 
reference to this specification and its drawings. It is contemplated that such 
modifications and alternatives are within the scope of this invention as subsequently 
25 claimed herein. 
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